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Emergent Calabi-Yau Geometry 
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We show how the smooth geometry of Calabi-Yau manifolds emerges from the thermodynamic 
limit of the statistical mechanical model of crystal melting defined in our previous paper In 
particular, the thermodynamic partition function of molten crystals is shown to be equal to the 
classical limit of the partition function of the topological string theory by relating the Ronkin 
function of the characteristic polynomial of the crystal melting model to the holomorphic 3-form on 
the corresponding Calabi-Yau manifold. 



Introduction: Stringy effects and quantum gravity ef- 
fects are expected to modify our concept of space and 
time, and understanding such effects quantitatively is an 
important step in the exploration of physics at the Planck 
scale. In the perturbative string theory, where the string 
coupling constant g s is assumed to be small, the string 
length is longer than the Planck length. In this regime, 
stringy corrections to spacetime geometry kicks in first 
and quantum gravity effects are obscured. In the context 
of the topological string theory on C 3 , it was pointed out 
in Ref. [2j that the partition function can be resummed, 
and the result can be expressed in terms of a statistical 
model of crystal melting for large g s . In our previous pa- 
per [H, we generalized this construction to an arbitrary 
toric Calabi-Yau 3-fold. In the crystal melting model, 
the Calabi-Yau geometry is discretized, and each atom 
of the crystal can be regarded as a fundamental unit 
of the geometry. In this Letter, we will show how the 
smooth Calabi-Yau geometry emerges from the discrete 
structure of the crystal melting model in the thermody- 
namic limit, where g s — > 0. The topological string theory 
is relevant for counting of microstates of black holes in 
the superstring theory [3j , and we expect that our result 
sheds some light on quantum nature of spacetime in the 
superstring theory also. 

A toric Calabi-Yau 3-fold M is a Kahlerjjuoticnt of 
C F+3 by U(l)® F , and its mirror manifold M is defined 
by the polynomial equation 0, [B} , 
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+ P(z,w)=0, (u,v)eC, (z,w)eC > 



Here P(z, w) is a Newton polynomial of the form, 



P(z,w) 



F+3 
i=l 



(1) 



(2) 



and Cj(£)'s are functions of the Kahler moduli t of the 
original toric 3-fold M. The exponents (nj,TOj) G Z 2 
correspond to lattice points of the toric diagram. For 
example, for the mirror of 0{— 1) + 0{— 1) bundle over 
CP 1 , P(z, w) is given by 



P(z, w) = 1 + z + w + e t zw. 



(3) 



In this Letter, we will show that the Newton polyno- 
mial P(z, w) for the mirror of a Calabi-Yau manifold is 
identical to the characteristic polynomial of the corre- 
sponding dimer model, which is the partition function of 
the model on a torus. The relation between P(z,w) and 
the characteristic polynomial had been discussed earlier 
in j6(. Here, we will prove their precise equality including 
the dependence on the moduli t. 

The dimer model enters into our discussion because 
of its relevance to counting of D brane bound states. 
Consider DO and D2 branes with a single D6 brane on 
M . The low energy effective theory is a supersymmetric 
quiver quantum mechanics characterized by a quiver di- 
agram on T 2 0,0]. In our previous paper [l| , we showed 
that the dimer model defined in Ref. [g, 0] gives a gener- 
ating function of the Witten indices for bound states in 
the quiver quantum mechanics [loj ]. 

We also constructed a statistical model of crystal melt- 
ing equivalent to the dimer model. The crystal consists 
of atoms of different types, each of which corresponds to 
a node of the quiver diagram for the quantum mechanics, 
and the chemical bonds between the atoms are dictated 
by edges of the diagram. The quiver diagram is drawn on 
T 2 . To construct the initial configuration of the crystal, 
we start with the universal covering of the diagram over 
the plane and pile up atoms on nodes following the rules 
prescribed in Ref. This initial configuration corre- 
sponds to a single D6 brane with no DO and D2 charges. 
Removing atoms from the crystal generates bound states 
with non-zero DO and D2 charges. Such molten crys- 
tal configurations are in one-to-one correspondence with 
perfect matchings of the dimer model defined in Ref. 0] , 
as shown in the Appendix of Ref. [l| . 

It is reasonable to expect that a classical geometric 
picture emerges in the limit of large DO and D2 charges 
since it can represent a large black hole in the super- 
string theory. The corresponding thermodynamic limit 
in the dimer model was studied in Ref. [ll|. In this 
Letter, we will show that the partition function of the 
dimer model evaluated in the thermodynamic limit is 
equal to the genus-0 limit of the partition function of 
the topological string theory on M. The dimer model 
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has been formulated to describe the non-commutative 
Donaldson-Thomas theory @, while the topological 
string theory for a general toric Calabi-Yau manifold is 
equivalent to the commutative Donaldson-Thomas the- 
ory [13], as shown by [lj|. A relation of our results to 
the wall-crossing formula 14, [ID, 17 1 between the 
non-commutative and commutative Donaldson- Thomas 
theories is currently under investigation. 

The emergence of Calabi-Yau geometry from the ther- 
modynamic limit has been observed in Ref. Q in the case 
of C 3 . In this Letter, we make the connection sharper 
and more explicit by showing the direct connection be- 
tween the partition functions of the crystal melting model 
and the topological theory for a general toric Calabi-Yau 
3-fold. 

Thermodynamic Limit of the Crystal Melting 
Model: The main object of study in this Letter is the 
partition function, 
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g s q -t qA 



(4) 



where fl(qo,qA) is the Witten index for bound states 
of qo DO branes and qA D2 branes on the A-th. 2-cycle 
(A = 1, dimi?2(M)) with a single D6 brane on a toric 
Calabi-Yau manifold. According to the dictionary in Ref. 
[l[, the Witten index is equal (up to a sign) to the num- 
ber of molten crystal configurations where qo is the total 
number of atoms removed, whereas the relative numbers 
of different types of atoms removed from the crystal are 
specified by qa's. Later we will identify g s as the topolog- 
ical string coupling constant and t as the Kahler moduli 
of the toric Calabi-Yau manifold M. 

The behavior of Z foro s — » can be evaluated by 



using the result of Ref. 11]. Consider a finite covering 



the original quiver diagram, N times in one direction 
and N times in another direction on T 2 . N is introduced 
as an infrared regulator, and we will take N — > oo at 
the end of the computation so that we have the dimcr 
model on the plane R 2 . The surface of the crystal is 
determined by the height function h over the plane. To 
define h, we start with the canonical perfect matching too 
of the dimer model corresponding to the initial crystal 
configuration with no DO or D2 charges. For any other 
perfect matching to, the superposition of Too and to gives 
a set of closed loops on the dimer graph. If to corresponds 
to a bound state with finite DO and D2 charges, to and too 
differ only in a finite region on the graph. The function hi 
(i: node of the the N x N cover of the quiver) is defined 
so that it is far away from the region where to and too 
differ, and it increases by 1 every time we cross a closed 
loop as we move inside of the region. The corresponding 
molten crystal configuration is obtained by removing hi 
atoms of the initial crystal over the node i. In particular, 



qo- 



(5) 



To take the thermodynamic limit, it is useful to intro- 
duce the Cartesian coordinates (x,y), < x,y < 1, on 
the N x N covering of T 2 . In the limit where g s <C 1 and 
1«JV, the height hi becomes a smooth function h(x, y). 
We rescale the height function by the factor of 1/N so 
that 



TV 2 / dxdy h(x, y) = 



qo 

N ' 



(6) 



to take into account the large N scaling of the partition 
function discussed in Ref. [181 ] and quoted as Theorem 2.1 
in Ref. 11]. The statistical weight in the thermodynamic 
limit is given by an integral of a surface tension a(dh), 
which is a function of the gradient of h, as (lij 



exp 



iV 2 maxf, / dxdy [—a(dh) — g s Nh(x, y)] 



(7) 



The integral of g s Nh(x,y) in the exponent comes from 
the weight factor e~ 3sqo in (0|, and we used In 
the thermodynamic limit, we look for a height function 
h(x,y) which maximize the exponent. 

To derive the macroscopic surface tension from the mi- 
croscopic crystal melting model, we first define a charac- 
teristic polynomial P(z,w) as the sum of perfect match- 
ings with weights assigned to edges of the dimer model 
on T 2 [ill ] . We then define its Ronkin function R(x,y) 
by 



R(x,y) 



lnP(e 



x+i6 e y+i<j>\ 



(2tt) 2 



(8) 



According to Theorem 3.6 of Ref. [ll|, the surface ten- 
sion a(d x h, d y h) is the Legendre transform of the Ronkin 
function with respect to (s, t) = (d x h,d y h) as 



R{x,y) 



-cr(s, t) + xs + yt] 



(9) 



The first step in relating the dimer model to the topo- 
logical string theory is to show that the characteristic 
polynomial P(z, w) of the dimer model is equal to the 
Newton polynomial P(z, w) for the mirror Calabi-Yau 
manifold (fT]), 



P(z,w) = P(z,w). 



(10) 



According to Ref. [20(, there is a one-to-one corre- 
spondence between perfect matchings of the dimer model 
on T 2 and bi-fundamental fields of gauged linear sigma 
model appearing in the Kahler quotient construction of 
the toric Calabi-Yau manifold M. They are then re- 
lated, by change of variables described in Ref. [2(|, to 
lattice points of the toric diagram and to terms z ni w mi 
in the Newton polynomial ^ . This shows that there is a 
one-to-one correspondence between terms in P(z, w) and 
P(z, w), as pointed out by Q. Furthermore we can show 
that their coefficients agree. The Kahler moduli of M are 
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the Fayet-Iliopoulos (FI) parameters of the quiver quan- 
tum mechanics. According to the dictionary between 
quiver gauge theories and dimer models 0], FI param- 
eters are associated with nodes of the quiver diagram, 
or equivalently to faces of the dimer model. Thus, we 
can identify the FI parameters with the magnetic fluxes 
through faces of the dimer model, which parametrize the 
energy of each perfect matching of the dimer model 11 1 . 
Each perfect matching appears in P(z, w) with the weight 
given by an exponential of the fluxes. On the other hand, 
the Newton polynomial is a sum of z ni w mi , each of which 
corresponds to a lattice point of the toric diagram and 
is weighted by an exponential of the Kahler moduli t Q • 
We have verified that the weight for perfect matchings 
and the weight for lattice points of the toric diagram 
agree, and this proves the identity (JTOj) . 

Combining ([7]) and ^ and discarding a term in total 
derivative in (x 7 y), which is justified by the subtraction 
of the linear piece in R(x, y) discussed in the next para- 
graph, we find 



exp 



iV 2 / dxdyR f ^-x, 



-y 



(ii) 



By rescaling (x,y) by the factor of gN/2, this becomes 



exp 



\ [ dxdyR(x,y) 
9 s J 



(12) 



Note that the N dependence has disappeared except that 
the range of the (x, y) integral has been rescaled by the 
factor of gN/2. For N — > oo with small but fixed g s , we 
have an integral over the whole (x,y) plane. 

The integral (fl"2"|) in the large N limit is divergent. To 
identify and subtract the divergent part, it is convenient 
to introduce the concept of the amoeba [2l| which is a 
subset of R 2 defined by, 



Amoeba = {(x,y) G R 2 : P{e x 



+i8 pV+H 







for some (6*, </>)}. 



(13) 



In the thermodynamic limit, the amoeba corresponds to 
the liquid phase of the crystal ll| . If there are no inte- 
rior points in the toric diagram, the complement of the 
amoeba is the solid phase, where the crystal retains its 
original shape. There, the Ronkin function R(x, y) is lin- 



ear 24]. If there are interior lattice points in the toric 
diagram, the amoeba acquires holes, inside of which are 
in the gas phase, where the Ronkin function is again lin- 
ear but the slope of the crystal surface is different from 
the original one. The integral (fT2|) becomes finite if we 
subtract the linear piece of the Ronkin function in the 
solid phase so that the partition function is normalized 
to be 1 for the initial crystal configuration. 

Topological String at Genus 0: Our next task is to 
compute the genus-0 topological string partition function 
Tq of the toric Calabi-Yau manifold M and compare it 



with the thermodynamic limit of the partition function 
of the crystal melting model (fT2j) . For this purpose, it 
is convenient to use the mirror Calabi-Yau manifold M 
defined by the equation ([1]) since can be evaluated by 
the classical period integral as, 



ho 



(14) 



where fi is the holomorphic 3-form on the mirror, and 
/3q is the Lagrangian 3-cycle which is the mirror of the 
6-cycle filling the entire toric Calabi-Yau manifold. 

According to the microscopic derivation of the mirror 
symmetry by Hori and Vafa [4j, the sigma- model on the 
toric Calabi-Yau manifold is equivalent to the Landau- 
Ginzburg model with the superpotential 

W(u,x,y) = e u P(e x ,e y ). (15) 

It was shown in Ref. [5J] that an integral of e w in a 3- 
dimensional subspace of the (it, x, y) plane can be trans- 
formed into a period integral of the holomorphic 3-form 
£1 on the mirror Calabi-Yau manifold. Thus, we should 
be able to evaluate (fl4|) as an integral of e w . To do so, 
we need to identify the contour of the integral. 

Since g s and t A, s in (01 are taken to be real in the 
dimer model, the Newton polynomial P(x, y) in our case 
is with real coefficients. The mirror manifold ([I]) has the 
complex conjugation involution, and thus the fixed point 
set is a natural candidate for In fact, following the 
mirror symmetry transformation as described in Ref. [4| , 
we find that the 6-cycle of the original toric Calabi-Yau 
manifold M corresponds to the real section in the (u, x, y) 
space in the mirror M. Thus, we find 



du 



dxdy e c 



'P(e 



'') - _ 



dxdy In P(e x ,e y ). 
(16) 

The divergent part of the integral in the {x, y)-plane can 
be removed by subtracting a linear term in lnP(e x ,e v ) 
for x,y — » oo as we did for (fT2"|) . The integral (fTT>)) is 
almost equal to the exponent of the partition function 
(fl"2)) of the crystal melting model, except that we do not 
have the averaging over the phases (8, <fi) to define the 
Ronkin function as in ((8|). It turns out that the integral 
over (x, y) in (|16|) removes the dependence on (9, </>), and 
thus the averaging process is not necessary. 

To see this, let us define a generalization of Tq for an 
integral of (x, y) with arbitrary phases as 



flM)-/**^,^). 



(17) 



Taking derivatives of T with respect to (0,(f>), we find 
the integrand becomes a total derivative in (x, y) as in 



= J dxdyi [a^-^p^j lnP(e x+te 7 e y+l4 '). 



(18) 
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If we choose (a, (3) so that it is not in the directions of the 
tentacles of the amoeba (|T3)) . the boundary term is re- 
moved by the regularization and we find (adg + (id^)T — 
0. Since a, (3 is arbitrary except in the directions of ten- 
tacles of amoeba, J- is independent of (9, </>) and agrees 
with its average. Namely, 

J-q = — J dxdybxP(x,y) = — J dxdy R(x,y). (19) 

Thus, we found that the thermodynamic limit of the 
partition function of the crystal melting model given by 
(fT2")l is equal to exp(— 4r-7"o), which is the genus-0 par- 
tition function of the topological string theory. This is 
what we wanted to show. 
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